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We have studied possible formation of vacua states in 
heavy ion collisions. Random phases of the chiral fields were 
evolved in a finite temperature potential, incorporating the 
breaking of Ua (1) symmetry. Initial random phases very 
quickly settle into oscillation around the values dictated by 
the potential. The simulation study indicate that an initial 
0=0 state do not evolve into aG^O state. However, an ini- 
tial 8/0 state, if formed in heavy ion collision, can survive, 
as a coherent superposition of a number of modes. 
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I. INTRODUCTION 

QCD in the chiral limit (m u = ma, = Tn s = 0) possesses 
a [^4 (1) symmetry. Spontaneous breaking of the symme- 
try requires a neutral pseudoscalar Goldstone boson with 
mass less than ^/%ra v , in addition to pion itself. However 
no such Goldstone boson is seen in nature. The problem 
was resolved by the discovery of non perturbative effects 
that violates the extra Ua{1) symmetry, 't Hooft showed 
that because of the instanton solutions of the Yang-Mills 
theory, the axial anomaly has nonzero physical effects, 
and there is not really a U(l) |j| . Term containing the so 
called vacuum angle (0) which breaks P and CP symme- 
try can be added to the Lagrangian. QCD requires very 
small 8 ~ 10~ 9 , which explain the apparent P and CP 
symmetries in strong interaction Dynamical breaking of 
the Ua(X) symmetry is also obtained in large N (color) 
limit of the SU(N) gauge theory (2|-|^]. In this approach 
the dominant fluctuations are not semi classical but of 
quantum nature. 

Recently Kharzeev, Pisarski and Tytgat |3 and 
Kharzeev and Pisarski |^| argued That, in heavy ion col- 
lisions, a non trivial vacua state may be created. In the 
limit of large number of colors, the axial Ua(1) symme- 
try of massless quarks may be restored at the deconfining 
phase transition. As the system rolls back to confining 
phase, it may settles into a metastable state with non- 
trivial 0. The idea is similar to the formation of disori- 
ented chiral condensate (DCC) Q. In DCC, space-time 
region is created, where the chiral condensate points in a 
direction, different from the ground state direction. Sim- 
ilarly, in a vacua states, space-time region with non- 
trivial 9 is created. 6 vacua state being P and CP odd, 
such a state, if produced in heavy ion collisions, may have 
some interesting signals, e.g. decay of rj into two pions, 



which is strongly forbidden in our world. 

After the suggestion of Kharzeev and co-workers HU , 
several authors have looked into various aspects of non- 
trivial vacua state, that may be formed in heavy ion 
collisions. p[-pd|. Buckley et al. JlC],[ll} numerically sim- 
ulated the formation of vacua states, using the effec- 
tive Lagrangian of Halperin and Zhitnitsky pi] ]. They 
assumed quenchlike scenario, rapid expansion of the fire- 
ball leave behind an effectively zero temperature region, 
interior of which is isolated from the true vacuum. Start- 
ing from an initial nonequilibrium state, they studied the 
evolution of phases of the chiral field, in a dissipative en- 
vironment. They saw formation of nonzero vacuum 
within a time scale of 10~ 23 sec. They choose to ig- 
nore the fluctuation-dissipation theorem, which require 
the dissipation to be associated with fluctuations. Re- 
cently we have performed a simulation study for vacua 
formation in heavy ion collisions The dissipative 

term was omitted. In a quench like scenario, random chi- 
ral phases were evolved in a zero temperature potential. 
It was shown that, for non-zero 0, randomly distributed 
initial phases oscillate about the values dictated by the 
potential. But, for the physically important case, the 
initial 0=0 state do not evolve into a9/0 state. 

Quench scenario is not expected to occur in heavy ion 
collisions. Unlike in a quench scenario, where the sys- 
tem is instantaneously brought to zero temperature, in 
heavy ion collisions, the fireball requires a finite time to 
cool. In the present paper, we have extended our study 
p^t to include the temperature dependence of the poten- 
tial. We have also added a dissipative term. Even in a 
quenchlike scenario, where the chiral phases are evolving, 
in essentially a zero temperature potential, a dissipative 
term is important. It can mimic the effects like interac- 
tion between the quarks, emission of pions etc. To be 
consistent with fluctuation-dissipation theorem, we have 
also included a fluctuation term (which was omitted by 
Buckley et al. |T^JTl[]). We represent it by a white noise 
source. In a realistic situation, if within a certain region, 
vacua state is formed, that region will be in contact 
with some environment (mostly pions etc.). Use of white 
noise source as the fluctuation term can take into account 
those interactions. In Rcf. |l2), we have considered a sin- 
gle event, which may be classified as a average type of 
event. vacua state being a rare event, it is not ex- 
pected to be produced in each and every event. In the 
present work, we have classified the events to single out 
the event where vacua state is most probable. As will 
be shown below, even in the most favored event, a initial 
= state, do not evolve into a6/0 state. 
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The paper is organized as follow: in section 1, we de- 
scribe the model. The results are discussed in section 2. 
The summary and conclusions are drawn in section 4. 



II. THE MODEL 

As in E2| , we have used the effective Lagrangian, devel- 
oped by Witten || , to study the 9 vacua state in heavy 
ion collisions. The effective non-linear sigma model La- 
grangian is ||, 
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£ = fn(^tr(d^d^U) + tr(M(U + U*)) 
-a(tr]nU - 9) 2 ) 



(1) 



where U is a 3 x 3 unitary matrix with expansion, 
U = U (l + iJ2 t a Tr a /U + 0(tt 2 )), Uo being the vacuum 
expectation value of U, t a are the generators of U(3), 
(Trt a t b = 5 ab ) and 7r a are the nonet of Goldstone boson 
fields. M is the quark mass matrix, which is positive, 
real and diagonal. We denote the diagonal entries as 
/if. They are the Goldstone boson squared masses, if 
the anomaly term a (a/N in ref. were absent. Be- 
cause M is diagonal U can be assumed to be diagonal , 

In terms of fa's, the potential is, 

Vfc) = E M? cos^i + a/2(J2 <f>i - ©) 2 ) (2) 

It may be noted that as ^2 (pi arose from tr lnf7, it is 
defined modulo 2tt. Vacuum expectation values of the 
angles </Vs can be obtained from the minimization con- 
dition, 



/i 2 sin < 



= o(0- 



(3) 



Witten || had discussed in detail, the solutions of 
these coupled non-linear equations. If the equation has 
only one solution then physics will be analytic as a func- 
tion of 9. The solution vary periodically in 9 with peri- 
odicity 2n H). Also, this solution must be CP conserving, 
whenever CP is a symmetry of the equation ||. However, 
it may happen that, EqJ^ has more than one solution. 
Then the solutions are not CP conserving, rather a CP 
transformation exchanges them. 

At finite temperature, the potential (Eq.||) is modified 
due to temperature dependence of the anomaly term a 
and also due to temperature dependence of the Goldstone 
Boson squared masses nf . In the mean field type of the- 
ory,temperature dependence of a and /i 2 can be written 



a(T) = a (l-T/T d ) 
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V 1 - T/Td 



(4) 
(5) 



where Td is the transition or the decoupling temperature. 
The subscript denote the values at zero temperature. 



In the present work we have used /j 2 = (114Mey) 2 , 
u3n = (161MeV0 2 ,M 2 o = (687) 2 and a = (492MeV) 2 
p,pl- With these parameters, the mass matrix in (0) can 
be diagonalised to obtain mj ~ 139 MeV, m v ~ 501 MeV 
and m,y ~ 983 MeV, close to their experimental values. 

Existence of metastable states (9 vacua states) can be 
argued as follows: the vacuum expectation values depend 
on the ratio a/fif. The ratio a/ [ij ~ (T d —T) 3 / 2 decreases 
with temperature. Then as the system rolls toward the 
chiral symmetry breaking, a/ [if may be small enough 
to support metastable states. In Q, it was shown that 
when a/ n\ < .2467 there is a metastable solution, which 
is unstable in ir° direction unless a < a cr , a cr / n\ ~ .2403. 

Appropriate coordinates for heavy ion scattering are 
the proper time (r) and rapidity (Y). Assuming boost 
invariance, equation of motion for the phases <^j's can be 
written as, 



d 2 ^ d 2 



^2/j 2 sin((t)i 



dx 2 dy 2 
-a(J2<fii-®))=<{T,x,y) (6) 



It is interesting to note that in this coordinate system, 
a dissipative term which decreases with (proper) time is 
effective. In addition, we have included the dissipative 
term with friction coefficient (rj) . rj take into account the 
interaction between the quarks. Expansion of the sys- 
tem, emission of pions etc. also contribute to r\. The 
exact form or value of r\ is not known. It is expected to 
be of the order of Aqc"£)=200 MeV. Temperature depen- 
dence may be assumed to be same as that of quarks and 
gluons. 77 for quarks and gluons generally shows a T 3 
type of temperature dependence Jl3]. We then assume 
the following form for 77, 



V (T)= VQ (T/T d f 



(7) 



with ?7o=200 MeV. In Eq.|, C is the noise term. We 
assume it to be Gaussian distributed, with zero average 
and correlation as demanded by fluctuation-dissipation 
theorem, 



< ( a {T,x,y) >= 
< Ca(r, x, y)Cb(r', x', y') >= 2^-5^ - r') 



(8) 



S(x - x')S{y - y')S ab (9) 

where a,b corresponds to phases <f> u , <f>d or <j) s . 

Solving Eq.^ require initial phases (0j) and their time 
derivatives (</>j) at the initial time (7^). We choose t\ 
=1 fm. <pi and (pi were chosen according to following 
prescription, 



(pi = R N / (1 + exp((r - r cr )/r) 
(pi = R N / (1 + exp((r - r cr )/T) 



(10) 
(11) 



where Rn is a random number within the interval 
[-27r/16,27r/16]. We use r cr =20 fm and T=.5 fm. 
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III. RESULTS 

The set of partial equations ([|) were solved on a 64 2 
lattice, with lattice spacing of a = 1/m and time inter- 
val of a/10 fm. We also use periodic boundary condi- 
tion. The cooling and expansion of the system was taken 
through the cooling law T(t) = T^/t, decoupling or 
transition temperature being Td =160 MeV. The cooling 
is rather fast and may not be suitable for heavy ion colli- 
sions. Still we choose to use it, to be as close as possible 
to a quench like scenario, which is ideal for vacua state 
formation. The phases were evolved for arbitrarily long 
time of 20 fm. We have performed simulation for 1000 
events. The idea is to see whether in any of these events, 
initial random chiral phases can evolve into a nontriv- 
ial O vacua state. We have considered two situation, (i) 
0=0 and (ii)0=47r/16, within the radius r cr and zero be- 
yond r cr . Physically, initial 0=0 is the most interesting 
case. In this case, emergence of non-zero fas will be the 
signature of non-trivial vacua state formation. It will 
show that during the roll down period, initially uncor- 
rected chiral phases get oriented in a certain direction. 
The 2nd situation, initial 0=4 7r/16, on the other hand, 
will explore the possibility of sustaining a initial vacua 
state in a heavy ion collision. 

A. Classification of Events 

vacua state, being a non-trivial state, should be a 
rare event. We do not expect it to be formed in each and 
every event. In the ensemble of 1000 events, vacua 
state may be formed in a few events only. It is then 
important to classify the events such that the rare events, 
where vacua state is formed, can be distinguished. For 
initial 0=0, emergence of non-zero fa is the signature 
of a nontrivial state. Then the space-time integrated 
value of the phases can be used to classify the events. 
Accordingly, in each event we calculate the classification 
parameter, 

G= [ Tdrdxdy \^fa\ (12) 

The event for which G is maximum, will be most prob- 
able event for forming a vacua state. In fig.l, we 
have shown the distribution of G for the 1000 events, 
for the two sets of calculations, (a) initial 0=0 and (b) 
initial 0=4-7r/16. Very interesting behavior is obtained. 
While for 0=0, the distribution falls with G, it is sharply 
peaked for 0=4 tt/16. We also note that for initial 
0=0, G ranges between [0.02-0.75]. It is apparent that 
the while the evolution depend on the initial conditions, 
in this case, the phases do not evolve to large values. 
Very small value of G, indicate that initial 0=0 state do 
not evolve into a 9 / state. In contrast, for initial 
0=47r/16, G ranges between [154.1-154.9]. The phases 



grow to large values in this case. Also very small differ- 
ence between the minimum and maximum values of G 
indicate that in this case, evolution of phases is largely 
independent of initial conditions. As will be discussed 
below, this is an indication that if formed, potential can 
support a vacua state. 

B. Evolution of spatially integrated phases 

In this section we study the temporal evolution of the 
spatially averaged chiral phases. For each trajectory, spa- 
tially averaged chiral phases, 

. , . I fa(x,y)dxdy . 

<4>i>=- 7T-— ,i = u,d,s (13) 

J dxdy 

were calculated. We also calculate the ensemble average 
of the spatially averaged phases, 

1 N 

< fa >en= ~yi<^ >j) i = U,d,8 (14) 

i=i 

In Fig. 2, for = 0, ensemble average of the spatially 
averaged phases < fa > en , along with its fluctuations for 
the 1000 events are shown (the black circles with error 
bars). All three phases, (f> u , fax and fa, show qualitatively 
similar behavior. For 0=0, zero temperature potential 
is minimized for <j> u — fai = fai—0. Even though fluctu- 
ations are large, ensemble averaged values remain zero 
throughout the evolution. The result clearly shows that 
on the average, initial 0=0 state do not evolve into a8/ 
state. In Fig. 2, solid lines shows the evolution of the 
spatially averaged phases (< fa >) in the most favored 
event. For initial 0=0, non-zero fa's will be indication 
of a vacua state. It is then expected that if a vacua 
state is formed, < fa >'s will remain nonzero for a sub- 
stantial time. We do not find such behavior. On the 
contrary, spatially averaged phases (< <f> u >,< fa% > and 

< fa, >) oscillate about zero, their vacuum expectation 
value. Amplitude of the oscillation decreases with time 
(due to friction). Even in the most favored event, initial 
= state do not evolve into a nontrivial state. 

In Fig. 3, same results are shown for initial 0=4 tt 
/16. The ensemble average of spatially averaged phases 

< fa >en, do not show large fluctuations. Very narrow 
range of G already indicated this. As a result, < fa > 
in the most favored event (the solid lines) closely fol- 
low < fa > en . For 0=47r/16, the minimization condi- 
tion gives, c/) u — 0.502, 0^=0.243 and <^> s =0.013, which are 
shown by the straight lines in Fig. 3, The spatially av- 
eraged phases are found to oscillate about these values 
throughout the evolution. As will be discussed below, 
continued oscillation of the phases indicate that in heavy 
ion collisions, an initial vacua state is sustained as a 
coherent superposition of a number of modes. 
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C. Correlation 



Wc define a correlation function, 



C(r) 



(15) 



such that the distance between the lattice points i and 
j is r. C(r) specifies how the three component 4> = 
{4>u, 4>d,4>s) at two lattice points are correlated. In a non- 
trivial O vacua state, phases will be correlated while in 
a trivial O vacua state (i.e. 0=0) phases will remain 
uncorrelated. In Fig. 4, we have shown the evolution of 
the correlation function for 0=0 and 0=47r/16 at differ- 
ent times, t=1,3,5,7 and 9 fm. The results are shown 
for the most favored event. Initially at r=l fm, for both 
the cases, correlation length is about 1 fm, the lattice 
spacing. Initial phases were random , uncorrelated. The 
correlation length increases at later times. But for 0=0, 
the increase is marginal. Thus initial uncorrelated phases 
do not develop correlation, confirming our earlier results 
that 0=0 state do not evolve into a non-trivial state. 
However for finite initial 0, correlation length increases 
very rapidly, and assumes quite large values. Physically, 
for finite 0, all the ^>'s try to align themselves in some di- 
rection (i.e. in the direction), thereby giving very large 
correlation length, even when a large distance separates 
them. This figure clearly demonstrated the possibility 
of parity odd bubbles formation in heavy ion collisions. 
Initially random phases evolve such that they points in 
the same direction, forming a large P, CP odd bubble. 



D. Momentum space distribution of phases 

Much insight to the process can be gained from the 
momentum space distribution of the phases. At each al- 
ternate time step, we apply a fast Fourier transform to 
the spatial data. The Fourier transformed data are then 
integrated over the angles to obtain momentum distribu- 
tion, 



1 f 2 * 
27r Jo 



, i = u,d,s 



(16) 



Each mode was averaged over some narrow bin. For 
initial = 0, evolution of ensemble averaged <jf>i(£;)'s at 
k=6.2 MeV (which we call zero mode) along with its 
fluctuations are shown (dots with error bars) in Fig. 5. 
Ensemble averaged </>i(fc)'s remain zero through out the 
evolution. They do not grow. The behavior is in agree- 
ment with the results obtained for ensemble average of 
spatially averaged (j) s 's. On the average, initial 0=0 state 
do not evolve into a ^ state. In Fig. 5, evolution of 
the zero mode (solid line) and higher modes k=18.5 and 
30.9 MeV (dash and dash-dashed line respectively), in 
the most favored event are also shown. The behavior of 



the modes in the most favored event is different. The 
zero mode as well as higher modes oscillates about the 
ensemble averaged zero value. Higher modes are largely 
suppressed compare to zero mode. The state is essential 
a zero mode state. The results again confirms that, even 
in the most favored event, a trivial state do not evolve 
into a non-trivial state. 

In Fig. 6, same results are shown for initial 0=47r/16. 
In contrast to 0=0 case, here the ensemble averaged 
zero mode quickly reaches some finite value, though it 
also shows large fluctuations. The modes in the most 
favored event show similar oscillatory behavior, oscillat- 
ing about the ensemble averaged value. And as before, 
higher modes are suppressed. However, suppression is 
not enough to label the state as a zero mode state. It 
is better to describe the state as a coherent superposi- 
tion of different modes. The results are in agreement 
with the evolution of spatially averaged 0j's. We found 
that spatially averaged </Vs execute oscillation about the 
value dictated by the potential. Present simulation re- 
sults are qualitatively different from the simulation re- 
sults of Buckley et al [jl0],[ll| • They found that for finite 
0, after a few oscillations, the zero mode quickly reaches 
the value dictated by the potential. Also higher modes 
are largely suppressed. The difference is essentially due 
to the strong dissipative environment (with out fluctua- 
tions), in which the phases were evolved. In the present 
case, friction assumes low value at late times, resulting 
in the continued oscillation of the phases, even at late 
times. 

Present simulations indicate that in heavy ion collision, 
initial random phases do not evolve into a non-trivial 
vacua state. However, an initial finite state, survive the 
evolution, as a coherent superposition of different modes. 
The result is encouraging, as it establishes that if a 
vacua state is created in a heavy ion collision, it will 
survive. This leaves open the possibility of creating a © 
vacua state in heavy ion collisions. Though our study 
indicate that initial random phases do not evolve into a 
non-trivial vacua state, it is not definitive. We have 
considered 1000 events only. Sample size may be too 
small to detect a vacua state, which, as noted earlier 
is a rare state. Also our event classification scheme may 
not be the best. It may be possible to devise a better 
criterion to distinguish between different events to detect 
a vacua state. 

What will be the signature of such a state. As such de- 
tecting states are difficult. Kharzeev and Pisarski |(| 
devised some observables to detect vacua state, and 
estimated the P-odd observables are of the order of 10~ 3 , 
a small effect. Also as discussed by Voloshin [[)) the so 
called signal of states may be faked by " conventional" 
effects such as anisotropic flow etc. Buckley et al pQ 11 
suggested the decay of 77(7/) — > -k + tt~ . It is strongly for- 
bidden in our world as C and P is violated. Thus in a 
P and CP odd bubble, the decay is a distinct possibil- 
ity. One may look for such decays, which are strongly 
forbidden in our world, but have definite probability in 
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a vacua state. Whatever be the signal of the vacua 
states, with a large number of modes contributing, sig- 
nal will be broadened, effectively diluting the detection 
probability. 
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IV. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the formation of 
vacua state in a heavy ion collision, in presence of dissi- 
pation and noise. Assuming boost-invariance, equation 
of motion of the chiral phases were solved on a 64 2 lat- 
tice, with lattice spacing of a=l fm. Simulation results 
for 1000 events were presented. A classification param- 
eter (space-time integrated phases) was used to distin- 
guish between different events. Evolution of spatially in- 
tegrated phases, different modes, indicate that , on the 
average, initial 0=0 state do not evolve into a finite © 
state. Even in the most favored event, a non-trivial 
state is not produced. The phase remain uncorrelated 
at later times also. However, as the sample space is not 
large, and as the vacua states are rare events, it is 
not possible to conclude definitely that vacua states 
will not be formed in heavy ion collisions. Situation is 
different if initially is nonzero. In this case, initial 
random phases quickly evolve to values dictated by the 
potential and continue to oscillate around it. Correlation 
studies also show that the initial random phases quickly 
develop large correlation, forming a large parity odd bub- 
ble. Fourier analysis of the modes indicate that initial 
6^0 states survive the evolution as a coherent super- 
position of a number of modes. 
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FIG. 1. Distribution of the classification parameter G, for 
6=0 and e=47r/16. 
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FIG. 2. Evolution of the ensemble averaged, spatially av- 
eraged chiral phases, cj> u , 4>d and <j> 3 for initial 0=0 (black 
circles with error bars). The solid lines show the evolution of 
the spatially averaged phases in the most favored event. 
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FIG. 3. Same as fig.2 for e=47r/16. 
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FIG. 4. Correlation function at different times for Q—0 and 
6=4 7r/16, in the most favored event. 




FIG. 5. Evolution of the ensemble averaged <j>i(k), k=6.2 
MeV (black circles with error bar), for initial = 0. The 
solid, dashed, dash-dashed lines are the evolution of (f>i(k) in 
the most favored event, for k~6.2, 18.5 and 30.9 MeV respec- 
tively. 




FIG. 6. Same as Fig. 5 for 9=4 tt/16 
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